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ABSTRACT 


Curved  wall  diffusers  designed  by  using  an  inverse  method  of 
solution  of  potential  flow  theory  have  been  shown  to  be  both  short 
and  highly  efficient.  These  features  make  this  type  of  diffuser 
attractive  in  thrust  ejector  applications.  In  ejectors,  however, 
the  flow  at  the  diffuser  inlet  is  nearly  a  uniform  shear  flow.  This 
paper  presents  a  method  used  in  examining  the  flow  velocity  along 
the  diffuser  wall  and  some  of  the  analytical  results  for  diffusers 
designed  with  potential  flow  theory  and  receiving  a  rotational  flow. 
The  inlet  flow  vorticity  and  the  diffuser  area  ratios  prescribed  in 
the  inverse  solution  of  the  irrotational  flow  are  the  parameters  of 
the  study.  The  geometry  of  a  sample  ejector  using  such  a  diffuser 
and  its  estimated  thrust  augmentation  ratio  are  also  presented. 
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INTRODUCTION 


The  existing  ’’Clemson  Inverse  Design  Program  for  Short  Curved  Wall  Dif¬ 
fusers"  ^  is  based  on  irrotational  flow  theory.  Much  experimental  data  have 
been  examined  for  inlet  flows  having  thin  boundary  layers  with  core  flows  well 
represented  by  potential  flow  theory.  In  such  cases,  experimentally  obtained 
wall  pressure  distributions  agree  with  the  theoretically  prescribed  distribu¬ 
tions  very  well.  In  ejectors,  however,  the  diffuser  inlet  flow  is  signifi¬ 
cantly  different  from  irrotational  flow. 

Based  on  the  inlet  velocity  measurements  reported  by  Hill  and  Gilbert^ 
it  appears  that  diffuser  inlet  flow  should  be  represented  by  a  shear  flow  of 
uniform  vorticity.  Under  this  condition  (admitting  a  shear  flow)  the  wall 
pressure  distribution  could  be  significantly  different  from  that  which  was 
prescribed  for  the  particular  diffuser  inlet  design.  Specifically,  this 
raises  a  concern  about  the  presence  of  an  adverse  pressure  gradient  or  a  de¬ 
celeration  in  diffuser  wall  velocity  which  may  result  from  the  inlet  shear 
flow.  A  moderate  deceleration  upstream  of  the  suction  slot  may  be  overcome  by 
increasing  the  suction  flow  rate.  A  strong  deceleration  will  result  in  flow 
separation.  A  combination  of  large  design  area  ratio  and  inlet  vorticity 
could  also  result  in  a  flow  reversal  at  the  diffuser  exit  because  of  the  shear 
flow  requirement,  even  in  an  inviscid  flow.  It  becomes  apparent,  therefore 
that  a  method  of  calculating  the  wall  velocity  of  a  shear  flow  in  this  type 
curved  diffuser  is  necessary  to  assure  no  adverse  pressure  gradient  along  the*  5 
diffuser  wall.  And  this  will  allow  us  to  obtain  the  high  thrust  augmentation^  ^ 
benefit  afforded  by  using  a  short,  curved  wall  diffuser  section  in  the  ejector*  p 
design.  ~  - 
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METHOD  OF  SOLUTION 


1.  Equations  of  Motion 

The  coordinate  system  is  defined  by  Figure  1. 
continuity  equation  is: 
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Define  the  stream  function  y(r,x)  such  that: 


For  axisymmetric  flow,  the 
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Then  equation  (1)  becomes: 
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In  addition,  V  satisfies  the  vorticity  equation  which  in  cylindrical  co¬ 
ordinates  for  an  axisymmetric  shear  flow  is: 
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Substituting  equations  (2a,b)  into  equation  (3)  yields  the  following 
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which  is  an  elliptic  partial  differential  equation. 


2.  Boundary  Conditions 

In  the  numerical  solution  of  equation  (4/  for  shear  flow,  we  used  the 
same  grid  network  as  for  irrotational  flow.  The  boundary  conditions  used  are 
as  follows. 


(1)  Inlet 

At  the  diffuser  inlet,  there  is  a  parallel  flow  in  the  x-direction, 
and  consequently  u  -  0*  The  y's  are  obtained  from  equations  (2a, b) 
as  follows 
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Therefore  y  only  changes  with  r,  hence 
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For  shear  flow, 
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Which  yields 


2  o,  in 


i*3 


From  equation  (5),  we  have: 
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where  R  -  r  <vwall) 


For  irrotational  flow  and  uniform  parallel  inlet  velocity 
have: 
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At  the  wall: 
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However,  fwali  may  not  necessarily  be  equal  to  (  ^irr )  wall  >  which 
means  that  for  the  same  diffuser  geometry,  flow  rates  for  the 
irrotational  and  rotational  (shear)  cases  may  not  be  the  same. 


Define  y  such  that 


and  then 


"iiAall  '  ’'’vail 


q.  c  Y(w  ,  -  t  BR.  ) 

nirr  o,in  3  m 


For  the  shear  flow  solution,  we  will  need  to  (i)  specify,  for  the 
inlet  boundary  condition,  the  values  of  wG  and  y  ;  (ii)  obtain  B 
from  equation  (8);  and  (iii)  specify  the  ’streamline  according  to 
equation  (5). 

Another  possibility  is  that  both  B  and  wQ  in  are  specified,  and  y  is 
to  be  calculated  from  equation  (8). 

In  the  inverse  design  program  if  we  have  pres*.  !bed  parallel 
but  non-uniform  irrotational  flows  at  the  inlet  and  the  exit,  then 
the  boundary  condition  at  the  inlet  would  be  as  follows: 


('"irrWl 


[l  'l  1 

|2  in  Wo , in  3 


copy  ’ 

msF&CTSBi 


SI  d 

**  j  2  -H 

c  t  c 

3.85*5 

m  *  t-*  o  •»  4 

Ji  f  ^ 
®  ;wo  c: 

5*  *—•*—*  r 
O  s-,  (_  „  .. 

I  I  65  «  D  ^ 


11-39 


r 


For  known  B  and  w  ,Y  can  be  calculated  from  the  following 
o,  in 
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and  the  streamlines  can  be  calculated  from  equation  (5)* 

(2)  Centerline 

Along  the  centerline,  we  have  ¥  =  0. 

(3)  Wall,  upstream  of  slot 

Along  the  wall,  upstream  of  the  slot,  we  have  f  =  Vuan  •  where 
^wall  same  as  the  one  calculated  at  the  inlet. 

(4)  Wall,  downstream  of  slot  and  inside  of  slot 

Let  3  and  3irr  be  the  fractions  of  the  flow  into  the  slot  for  the 
rotational  and  irrotational  cases  respectively.  These  3  1  s  are  re¬ 
lated  to  the  stream  functions  as  shown  in  the  following  expres¬ 
sions. 
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3  is  specified  as  an  input  to  the  analysis.  The  y  values  for  lines 
AB  and  AC  in  Figure  2  are  determined  from  equation  ( 10) • 

(5)  Diffuser  Exit 

Letting  Qirr  *  Volumetric  flow  rate  at  the  diffuser  exit  for  irrota¬ 
tional  flow  and 


Q  *  Volumetric  flow  rate  at  diffuser  exit  for  shear  flow 


results  in: 
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Here  B  is  specified,  f .  and  R are  obtained  from  the  Clerason 
Inverse  Solution  Computer  Program,  and  wq  is  the  center  line 
velocity  at  the  diffuser  exit.  (R  MCTArTT,e  is  the  radius  at  the 
stagnation  point  A  of  Figure  2.) 

(6)  Exit  of  Slot 

At  the  slot  exit,  a  shear  flow  having  vorticity  B  is  assumed. 


Referring  to  Figure  4: 


or 
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The  following  development  is  to  determine  the  stream  function  in  the 
slot. 
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But  from  Figure  3  we  see  that,  q  -  q q  ^  -  By 
and  therefore 
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The  yfs  at  the  exit  of  the  slot  can  be  prescribed  using  equation 

(R-R  ) 
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3.  Computational  Procedures 
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(1)  For  irrotational  flow,  prescribe  9  irr  parallel  at  inlet  and  exit  of 
diffuser.  Prescribe  a  slightly  increasing  velocity  distribution  at 
wall,  upstream  and  downstream  of  the  slot,  with  desired  deceleration 
across  the  slot  region. 

(2)  Solve  the  inverse  problem  for  irrotational  flow.  The  grid  work  in 
terms  of  x*s  and  rfs  and  the  velocity  distributions  at  inlet  and 
exit  of  the  diffuser  and  at  the  slot  exit  are  then  obtained. 

(3)  Solve  the  rotational  flow  equation  for  V  using  the  previously  dis¬ 
cussed  boundary  conditions. 

(4)  Solve  for  the  velocity ' distribution  at  the  wall.  If  the  velocity 
distribution  at  the  wall  is  slightly  increasing  with  an  abrupt 
deceleration  across  the  slot,  the  diffuser  geometry  is  satisfactory. 
If  there  is  a  deceleration  along  the  wall,  change  the  irrotational 
wall  velocity  distribution  and  go  back  to  (2). 
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Using  equations  (2a,b)  we  may  compute  the  velocity  distribution  along  the 
wall  from  the  t  distribution*  i.e.,  with 

- 

In  the  wall  velocity  computation,  the  vorticity  is  to  be  specified  as  one 
of  the  parameters  at  the  diffuser  inlet,  and  in  ejector  application  this 
parameter  is  determined  from  the  mass  ratio  of  the  primary  and  the  secondary 
flow  rates  and  the  geometry  of  the  mixing  chamber.  Using  this  method, 
computations  were  carried  out  on  the  University's  IBM  370/3033  digital 
computer.  Before  a  more  generalized  inverse  design  program  is  formulated,  the 
method  outlined  above  will  be  used  to  examine  the  wall  velocity  distribution 
to  assure  the  absence  of  any  adverse  pressure  gradient  in  ejector  design.  A 
more  detailed  version  of  the  analysis  will  be  published  as  Mr.  Francois 
Ntone's  thesis  for  a  Master  of  Science  Degree  in  Mechanical  Engineering  at 
Clemson  University,  which  is  expected  to  be  completed  by  May  1981.  Results  of 
the  sample  computations  using  the  above  outlined  analysis  are  presented  in  the 
next  section. 


ANALYTICAL  RESULTS 


Figure  5  represents  an  axisyrametric  diffuser  having  an  area  ratio  of  2.5 
to  1.  Figure  6  shows  the  velocity  distributions  along  the  diffuser  walls. 
The  circles  denote  the  velocity  distribution  prescribed  to  the  "Inverse  Design 
Program"  for  irrotational  flow,  and  the  triangles  represent  the  velocity  dis¬ 
tribution  along  the  diffusers  depicted  in  Figure  5*  The  second  velocity  dis¬ 
tributions  were  computed  for  inlet  shear  flow  having  a  nondimensional  vortic¬ 
ity  value  of  0.33  where  velocity  and  length  are  normalized  by  inlet  center 
line  velocity  and  inlet  width  respectively.  It  is  apparent  that  even  though 
the  prescribed  diffuser  wall  velocity  distribution  upstream  of  the  slot  has  a 
strong  acceleration,  the  wall  velocity  distributions  may  experience  deceler¬ 
ation  in  the  diffuser  when  a  shear  flow  is  admitted.  In  actual  operation  a 
flow  separation  would  most  likely  take  place. 

Figure  7  shows  three  velocity  distributions  for  the  diffuser  depicted  in 
Figure  5;  these  are  specifically  for  inlet  vorticity  values  of  0.33,  0.5,  and 
0.65.  There  is  a  significant  change  in  the  magnitude  of  the  velocities,  but 
their  gradients  vary  only  slightly.  As  expected,  larger  vorticity  results  in 
more  deceleration.  The  larger  change  of  wall  velocity  magnitude  results  from 
the  assumption  in  our  analysis  that  the  diffuser  center  line  velocity  is  kept 
constant,  therefore  the  wall  velocity  decreases  more  when  the  inlet  shear  flow 
vorticity  increases  more.  Judging  from  the  distributions  of  this  figure,  a 
flow  separation  would  likely  take  place  in  all  three  cases.  Figure  8  shows 
the  wall  velocity  distributions  for  the  diffuser  depicted  in  Figure  5,  under 
the  conditions  of  (i)  admitting  an  inlet  shear  flow  with  a  vorticity  of  0.5 
and  (ii)  three  different  slot  suction  rates.  A  suction  rate  Of  5.6%  was  the 
design  value.  At  a  suction  rate  of  8.5%,  the  velocity  gradient  upstream  of 
the  slot  has  not  improved  significantly  in  comparison  to  the  distribution  re¬ 
sulting  from  the  5.6%  suction  flow.  Perhaps  none  of  these  distributions  will 
yield  a  fully  attached  flow  throughout  the  diffuser. 


Figure  9  shows  a  diffuser  designed  for  irrotational  flow  and  having  an 
area  ratio  of  2.5  to  1  with  its  suction  slot  slightly  modified.  Figure  10 
shows  the  rotational  velocity  distribution  along  the  walls  of  diffusers 
depicted  in  Figures  1  and  5  with  inlet  vorticity  of  0.5.  Apparently,  there  is 
no  significant  difference  in  wall  velocity  distribution  for  these  two  shear 
flow  inlet  cases.  Figure  11  shows  a  diffuser  designed  for  irrotational  flow 
and  having  an  area  ratio  of  2  to  1.  Figure  12  shows  wall  velocity  distribu¬ 
tions  for  the  diffuser  depicted  in  Figure  11.  The  velocity  distribution 
represented  by  circles  is  for  irrotational  flow,  and  that  represented  by  tri¬ 
angles  is  for  shear  flow  with  velocity  of  0.5.  In  this  case  hardly  any  decel¬ 
eration  is  detected  for  the  rotational  flow,  therefore,  one  should  expect  a 
highly  effective  diffuser  even  when  a  shear  flow  inlet  condition  is  imposed  at 
the  curved  wall  diffuser  inlet.  Figures  13  and  14  are  similar  to  Figure  12 
and  are  for  an  area  ratio  of  1.5  to  1. 


SAMPLE  EJECTOR 


Now  the  variation  of  wall  velocity  distribution  due  to  vorticity  is 
understood,  and  one  may  select  ejector  design  parameters  in  such  a  way  that  no 
deceleration  takes  place  along  the  diffuser  wall  and  the  diffuser  flow  can  be 
maintained  attached  without  excessive  boundary  layer  «suction. 

Figure  15  shows  an  example  of  an  air-to-air  ejector  for  thrust  augmenta¬ 
tion  with  the  area  ratio  of  the  secondary  flow  to  the  primary  flow  at  the  mix¬ 
ing  chamber  inlet  of  40  to  1,  and  the  mass  ratio  of  the  secondary  flow  to  pri¬ 
mary  flow  of  12  to  1.  This  ejector  has  a  static  pressure  of  -1.25/psig  at  the 
exit  of  a  mixing  chamber.  From  this  pressure  level  one  may  use  a  highly 
effective  short  diffuser  with  an  area  ratio  of  2.2  to  achieve  an  exit  pressure 
of  the  ambient  level. 

Based  on  the  conventional  definition  of  thrust  augmentation  ratio  <p  , 
namely;  the  momentum  of  the  mixture  of  the  primary  and  secondary  flows  at  the 
ejector  exit  to  the  momentum  of  the  primary  flow  at  the  mixing  chamber  inlet, 
this  sample  ejector  has  a  vlaue  of  <p  of  2.182.  Since  the  short  diffuser  of 
the  sample  ejector  uses  boundary  layer  control,  an  auxilliary  ejector  will  be 
used  to  achieve  the  necessary  boundary  layer  suction.  The  removed  air  will  be 
discharged  in  the  same  direction  as  the  thrust  ejector  and  thus  contributes  to 
the  overall  thrust.  Therefore  a  modified  thrust  augmentation  ratio  <f>2  »  de¬ 
fined  as  the  total  momentum  at  the  ejector  exit  to  the  sum  of  the  momentum  of 
the  primary  flows  at  the  mixing  chamber  inlet  and  the  suction  chamber,  must  be 
used  for  this  type  of  ejector.  We  found  that  for  the  sample  ejector,  $2  has  a 
value  of  1.830  provided  that  a  suction  fraction  of  10%  is  adequate  for  the 
boundary  layer  control,  and  the  mass  ratio  of  the  auxilliary  ejector  of  4  to  1 
can  be  accomplished.  These  assumptions  are  considered  to  be  reasonable.  It 
is  noteworthy  that  the  overall  length  of  the  ejector  is  much  reduced  relative 
to  the  ones  having  diffusers  designed  by  using  the  concept  of  incipient  sepa¬ 
ration.  The  performance  level  in  thrust  augmentation  ratio  compared  favorably 
to  those  reported  in  reference  3. 
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Figure  1.  Coordinate  System  and 
Velocity  Components 


Figure  2.  Wall  Streamlines  and 
the  Branching  Point 


Figure  3.  Velocity  Distributions  for  Shear  Flow 


Figure  4.  Coordinate  System  at  Slot  Exit 
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Figure  5.  Computer  Plot  of  <f>  -  t  Grid  for  an  Axisymmetric  Diffuser,  AR 
Suction  Rate  *  5.6% 
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Figure  6.  Velocity  Distributions  for  Rotational  and  Irrotational  Flow 
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Figure  7.  Velocity  Distributions  for  Three  Rotational  Flows  with  Suction  Rate 
of  5.6% 
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Figure  8.  Velocity  Distributions  for  Rotational  Flow  with  B  -  0. 5  and  Three 
Suction  Rates 
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Figure  14,  Velocity  Distributions  for  Rotational  and  Irrotational  Flows,  B  ■ 
0.5,  Suction  Rate  -  8.5% 


ample  Air-to-Air  Ejector  with  10%  Contraction  for  Mixing  Chamber 
Deceleration  Along  Diffuser  Wall 


